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We study magnetic phases of two-component mixtures of ultracold fermions with repulsive inter- 
actions in optical lattices in the presence of hopping imbalance. Our analysis is based on dynamical 
mean-field theory (DMFT) and its real-space generalization at finite temperature. We study the 
temperature dependence of the transition into the ordered state as a function of the interaction 
strength and the imbalance parameter in two and three spatial dimensions. We show that below the 
critical temperature for Neel order mass-imbalanced mixtures also exhibit a charge-density wave, 
which provides a directly observable signature of the ordered state. For the trapped system, we 
compare our results obtained by real-space DMFT to a local-density approximation. We calculate 
the entropy for a wide range of parameters and identify regions, in which mass-imbalanced mixtures 
could have clear advantages over balanced ones for the purpose of obtaining and detecting quantum 
magnetism. 
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Ultracold atoms in optical lattices have become in 
recent years very convenient and powerful systems for 
quantum simulations of condensed matter materials [l|- 
y]. In a variety of cases, these gases show even richer 
physics than soHd state systems. While spectacular 
progress has been made in realizing different many-body 
phenomena in the past few years, many quantum phases 
present in condensed matter physics still represent very 
challenging goals for cold gas experiments. One of these 
challenges is to obtain quantum magnetism in optical lat- 
tices. Considerable progress towards this goal was made 
in a recent experiment [4] with ultracold ^^K prepared in 
a mixture of two different hyperfine states. 

Here we focus on two-component mixtures of repul- 
sively interacting fermions with an imbalance in the 
hopping amplitudes that can be realized, for example, 
by loading ultracold ^Li and ^^K [5] or different alka- 
line earth atoms [6] into an optical lattice, or by using 
spin-dependent optical lattices for two different hyperfine 
states [7]. We consider a Fermi-Hubbard Hamiltonian of 
the following type: 

n =-tA E(«l% + H-c.) - tB Y^i^lbj + H.c.) 

-^U^hiATliB^^ ^ {Vi - IJ.a)nia, (1) 

i i a=A,B 

where Ia and ts are the hopping amplitudes of fermionic 
species A and 5, respectively. The notation (z, j) indi- 
cates a summation over nearest-neighbor sites, and U 
is the amplitude of the on-site repulsive {U > 0) in- 
teraction of different species with corresponding densi- 
ties fiiA and fiiB' In the last term, Vi is the ampli- 
tude of the external (e.g., harmonic) potential at lat- 
tice site z, and jUa is the chemical potential of species 
a (henceforth we restrict ourselves to the case jj.a = 
Mb = m)- The Hamiltonian ([1]) corresponds to the 



single-band approximation; in other words, we consider 
the case of a sufficiently strong lattice potential, V\at ^ 
5Er. Thus, the hopping amplitudes ta are given by 

[8] , where 



t 
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2V^i:Vi?. 



Era = h'^k'^/2ma is the recoil energy, k is the wave num- 
ber determined by the wavelength of the laser forming 
the optical lattice, and m^ is the mass of species a. The 
amplitude V^[^^ of the lattice potential can be different for 
the two components V^l^^ ^ ^ 1^^^ ^ , which results in the 
possibility to realize an imbalance in the hopping am- 
plitude even for different hyperfine states of same atom 
{rriA = tub) [7]. 

By using the Schrieffer-Wolff transformation in the 
limit tA,B "C U near half filling, nAi-\-TiBi ^ 1, the Hamil- 
tonian ffl can be mapped onto an effective spin Hamil- 
tonian [qI, 110| . For the system under study, this transfor- 
mation results in the anisotropic Heisenberg model. 



^eff = J\ 









^5f 



SrsJ), (2) 



with coupling constants J± = 4:tAtB/U and 
J\\ = 2{t\^t%)/U. Note that here and below, all 
"magnetic" characteristics refer to a pseudospin made 
of two different species {A and 5), which corresponds 
to the spin operators S^ = {alai — h\hi)/2 and 
{S^ + iSj) = dlbi. Thus it should be mentioned that 
the anisotropy in the spin model is not related to the 
actual spatial directions or quantization axes in the 
original optical lattice. 

In the presence of hopping imbalance, i.e., tA i^^B^ the 
antiferromagnetic coupling Jy (in Z direction) is always 
larger than the coupling Jx (in the XY plane). Note 
that in the limit of large hopping imbalance, \tA—'^B\ -^ 
{tA + ^b), we find high anisotropy, Jy :^ J_\_; thus, the 
second term in the Hamiltonian (|2]) can be considered 



as a perturbation, and we arrive at the Ising model. 
Therefore, the imbalance in the hopping amplitude can 
be treated as a parameter that effectively reduces quan- 
tum fluctuations of the system under study. As shown 
below, this reduction results, in particular, in a relative 
increase of the Neel-ordered region even for larger tun- 
neling ta ~ U. 

Note that the above mapping to the effective spin 
model is presented in order to give a better physical un- 
derstanding of the system under consideration. Natu- 
rally, for quantitative theoretical predictions in the in- 
termediate coupling region {ta ~ U) it is necessary to 
use nonperturbative numerical methods, in combination 
with the full Hamiltonian ([it. B elow we apply dynamical 
mean-field theory (DMFT) [11| to examine the influence 
of hopping imbalance on the ordered phases in mixtures 
of ultracold repulsive fermions. 

Before diving into details, let us introduce parameters 
that are useful for further analysis: the average hop- 
ping t = {tA + ^b)/2 and the dimensionless hopping- 
imbalance parameter At = {tA — tB)/{tA + ^b)- Unless 
otherwise specified, we consider atoms of type A to be 
lighter than atoms of type B {tA > ^b), such that the 
hopping-imbalance parameter is positive. At G [0, 1]. 

To study bipartite order, such as antiferromagnetism, 
it is necessary to consider a two-sublattice configuration 
within the self-consistency conditions of DMFT in the 
homogeneous case [12|. Note that physical observables 
corresponding to these sublattices, in contrast to the 
case of a balanced mixture {tA = ts = t), cannot be 
directly associated with quantities corresponding to the 
two different species. Instead, the Green's functions are 
defined as Gsa{i^n) = Csa j D{e)/{CiaC2a - ^^)de [lH, 
where s = {1,2} and s are the sublattice index and its 
opposite, respectively, C,sa = ^^n + Ma — ^sa{'i^^n)i the 
Matsubara frequency corresponding to the temperature 
T (we set ks = 1) is ujn = 7r(2n + 1)T with n G Z, the 
self-energy T^sa{i^n) carries information about local in- 
teractions, and D{e) is the density of states [we present 
results for D{€) corresponding to cubic or square lattice 
geometries below]. To solve the impurity problem, we 
use the exact diagonalization solver [11|. We also verify 
our results for selected parameters by performing calcu- 
lations with other highly accurate numerical algorithms, 
such as the continuous- time Monte Carlo [13| and numer- 
ical renormalization group [14] impurity solvers. 

For the homogeneous system, we concentrate on the 
case of half filling, fi = U/2. At low temperatures its 
ground state exhibits Neel (antiferromagnetic) ordering. 
The corresponding phase diagram obtained by DMFT 
for the cubic lattice is depicted in Fig. [H One can 
conclude that the hopping imbalance results in a rel- 
ative increase of the critical temperature (note that it 
is scaled in units of average hopping t). In the atomic 
limit {U ^ t), one can derive this from the mean- field 
value of the Neel temperature in the Heisenberg model. 
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FIG. 1. (a) Critical temperature for Neel ordering of the half 
filled Hubbard model at different values of the mass- imbalance 
parameter in a cubic lattice, obtained within DMFT. (b), 
(c) Contour plots indicating the absolute value of the order 
parameter mstag = \nAi — nBi\ for the balanced (b), At = 0, 
and maximally- imbalanced (c), At = 1.0, mixtures. 



Tn = 6JS{S ^l)/3 [15|, where S is the fermions' spin. 
Taking J = J\\ = 2(t^ + ^b)/^' ^^^ obtains for constant 
U/t the relation Tat (At) /Tat (0) = 1 + At^ 

We also perform analogous calculations for a two- 
dimensional system with square lattice geometry. We 
would like to emphasize here that the Mermin- Wagner 
theorem, which excludes long-range order at finite tem- 
perature in low dimensions {d < 2) for continuous sym- 
metries, is applicable only to the rotation symmetry of 
the Hamiltonian (J2j) in the XY plane. Thus, for Jy > J^, 
long range Z-antiferromagnetic order is permitted [15]. 
One can verify that, even for d = 1 [16], the excitation 
(magnon) spectrum has an energy gap A oc (Jy — J±). 
This protects the system from the infinite number of low- 
energy excitations, which destroy antiferromagnetic or- 
der in the balanced (t^ = t^) system at T > in d < 2. 
Regardless of whether hopping imbalance is present or 
not, a critical temperature in d = 2 can be defined if one 
uses a mean- field approach (or DMFT, as applied here). 
Our results for the dependence of the transition tempera- 
ture on the imbalance parameter and interaction strength 
are similar to the dependence for that of the cubic lattice 
as shown in Fig. [1] but with lower values of the critical 
temperature; e.g., we obtain max[T^^^(At = 0)] ^ 0.39t 
and max[T2^(At = 1)] ^ 0.45t. 

It should be noted that the data obtained for the tran- 
sition temperatures in the limiting case of a mixture with 
equal hopping amplitudes (At = 0) quantitatively agree 
with results presented in Refs. [13, llll, where DMFT 
was also used. We emphasize that it has been well- 
established that, in the balanced case, DMFT overesti- 
mates the critical temperature for moderate to large in- 
teractions {U/t > 5) in comparison with exact quantum 
Monte Carlo (QMC) simulations [19|,|20| and a dynamical 
cluster approximation (for quantitative comparison see 
[2l|). However, in the limit of large interactions, existing 
calculations [22] show that the qualitative increase of Tc 
due to an imbalanced mixture will in fact be enhanced 
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FIG. 2. Phase diagram of a two- component fermionic mixture 
obtained by DMFT at finite temperature. The temperature 
is scaled proportionally to the average hopping, T = 0.4t. 
The contour plots to the right indicate the absolute value of 
the staggered magnetization ?7istag, the charge- density wave 
parameter c — \D — K\l {D ^ K), and the average double oc- 
cupancy D, where the average is taken over the sublattices 
s = 1,2. 



by the use of an exact method such as QMC calcula- 
tions. We expect that this is also true for intermediate 
interactions. 

One can also gain important information from a 
phase diagram, where the hopping amplitudes of dif- 
ferent species are scaled independently (for bosons, see 
Refs. [23, '24^), which is presented in Fig. [21 It contains in- 
formation on both magnetic ordering and the conductive 
(localization) properties of the system. The Neel phase 
is characterized by an order parameter, the staggered 
magnetization, equivalent to a spin-density wave. How- 
ever, from Fig. [2] it is clear that in the mass-imbalanced 
mixture one has an additional type of order, a charge- 
density wave (CDW), which, in general, is nonzero in- 
side the ordered phase except for equal hopping ampli- 
tudes, tA = ts (balanced mixture), where it vanishes. 
To define whether the system is in the metallic or insu- 
lating state at half filling, we also analyze both the dou- 
ble occupancy, Di = {fiAiTiBi), and the zero occupancy, 

Ki = {{l-hAi){l-nB^)). 

The origin of the charge-density wave can be more 
clearly understood from the analysis of spatial distribu- 
tions. To this end, we apply an extension of DMFT, 
the real-space dynamical mean-field theory (R-DMFT) 
[25|, |26| . It effectively captures long-range order together 
with effects originating from a finite size of the sys- 
tem and an external trapping (harmonic) potential. In 
Fig. Owe show real-space distributions of the total filling 
^i = nAi-\-TiBi and the occupation difference (magnetiza- 
tion) rrii = riAi — tlbi both in two- and three-dimensional 
optical lattice geometries. 

From Figs. [3fb)-(d) we see that the emergence of a 
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FIG. 3. Upper row (a)-(d): Magnetization and total den- 
sity distribution of mass-imbalanced repulsive fermions in a 
harmonic trap at different temperatures (R-DMFT, square 
lattice, d — 2). (e) The same distributions obtained within 
the local-density approximation (LDA+DMFT, lines) and R- 
DMFT (dots, corresponding to individual sites) at tempera- 
ture T = 0.2tA for a cubic lattice in d = 3. Other param- 
eters used in (a)-(e): ts = 0.5tA, U = lOtA, M = t//2, and 
Vi — 0.1tA(n/a)^, where a is the lattice constant. 



charge-density wave at half filling {ui ^ 1) corresponds to 
the fact that in the Neel-ordered state the sites z, which 
are occupied by the heavier component, have a larger 
total filling (n^ > 1); i.e., they have an enhanced double 
occupancy Di {Di > Ki). At the same time, the sites j 
that are mainly occupied by the lighter component have a 
smaller total filling (rij < 1); i.e., they have an enhanced 
zero occupancy Kj {Kj > Dj) due to a higher mobility of 
the lighter atoms. One could possibly use this property as 
an alternative signature of magnetic ordering of ultracold 
hopping-imbalanced fermions in optical lattices, e.g., by 
performing site-resolved detection of atom pairs. 

It should be noted that in the real-space density dis- 
tributions of mass-imbalanced mixtures with fiA = Mb 
there is always a "ferromagnetic" ring present [it is di- 
rectly visible in Fig. [3l^a)]. This ring originates from the 
larger kinetic energy of the lighter component that, in 
turn, results in a wider spatial distribution of these atoms 
in a trap. We conclude from Figs.[3l^b)-(d), that this ring 
remains present at temperatures above the critical tem- 
perature and, thus, cannot be uniquely identified with 
quantum magnetism. 

In Fig. [3l^e), we present a comparison of two computa- 
tional approaches: R-DMFT and DMFT combined with 
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FIG. 4. Phase diagrams of the half filled Hubbard model 
(/xo = ^/2) for balanced (a) and imbalanced (b) mixtures 
obtained within DMFT in the homogeneous (nontrapped) 
case. Shaded areas correspond to the Neel (antiferromag- 
netic) phase. The colored and dark- gray lines are isentropic 



a local-density approximation (LDA). We find very good 
agreement in the center of the trap (including the mag- 
nitude of the charge-density wave) and at the edges. But 
in the intermediate region between the antiferromagnetic 
core and the "ferromagnetic" ring, LDA fails to repro- 
duce the detailed structure as it does not account for 
a proximity effect; i.e., it neglects the influence of sur- 
rounding sites with different densities. In contrast, R- 
DMFT shows a wider spatial range of stability of the 
antiferromagnetic order that results in an interesting in- 
terplay between antiferromagnetic correlations and the 
"ferromagnetic" ring. Note that proximity-induced an- 
tiferromagnetic order in a trapped system has also been 
shown previously in R-DMFT studies of balanced mix- 
tures [25|,|26[. 

Finally, as our approach allows for the extraction of 
expectation values of the local filling with high precision, 
one can also perform an entropy analysis for the homoge- 
neous system. Using a Maxwell relation for the entropy s 
per lattice site i, ds/dfi = dn/dT^ where n = ^^ nf^ /2^ 
one obtains s{iiq^T) = J^° {dn/dT)dfi. Therefore, by 
performing calculations at different values of the interac- 
tion strength U and temperature T, we obtain the isen- 
tropic curves presented in Fig. SI 

As one can conclude from Fig. IH the approach towards 
quantum magnetism in optical lattices based on prepa- 
ration of the atomic mixture in the Fermi-liquid state, 
and consequent adiabatic increase of the interaction also 
works for the case of mass-imbalanced fermions. The re- 
sults obtained from our DMFT calculations show that, 
in order to cross the transition boundary, one needs to 
prepare the system in a state with entropy per particle 
s < In 2 [27], independently of the value of the mass- 
imbalance parameter. However, if one starts with higher 
entropy values, then it is more favorable to use imbal- 



anced mixtures as they allow a much closer approach to 
the critical region by adiabatic increase of the interac- 
tion. For example, at 5 = 0.75 one obtains Tmin = 0.82t 
for At = 0, while at the same entropy value Tmin = O.blt 
for At = 0.8. One can also see that at the same rela- 
tive values of the temperature in the Fermi-liquid region 
{U ^t) the entropy increases with hopping imbalance. 

It should be noted that Fig. 2] is obtained for the ho- 
mogeneous (i.e., nontrapped) case. Thus, the depicted 
isentropic curves include only one of the possible mecha- 
nisms of cooling, namely the Pomeranchuk effect in two- 
component ultracold fermionic mixtures [17[ . Let us em- 
phasize that, in addition to this effect, the trapped sys- 
tem can also be cooled by a redistribution of atoms within 
the trap (see Refs. [28-30| for details). 

To summarize, in this Letter we demonstrate advan- 
tages of hopping-imbalanced two-component mixtures of 
ultracold fermions for approaching quantum magnetism. 
We compare critical temperatures and isentropic curves 
obtained by DMFT and show that for imbalanced mix- 
tures the critical region can be approached at higher en- 
tropy values than for balanced mixtures. We also ob- 
serve that hopping-imbalanced ultracold Fermi gases in 
the Neel-ordered state have additional long-range order 
in the form of a charge-density wave, which provides an 
alternative signature for detecting quantum magnetic or- 
der in optical lattices. 
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